Chiral spin states in the pyrochlore Heisenberg magnet 
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Fermionic mean-field theory and variational Monte Carlo calculations are employed to shed light 
on the possible uniform ground states of the Heisenberg model on the pyrochlore lattice. Among 
the various flux configurations, we find the chiral spin states carrying ±7r/2 flux through each 
triangular face to be the most stable, both within the mean-field theory and the projected wave 
function studies. Properties of the spin-spin correlation function and the chirality order parameter 
are calculated for the projected wave functions. Mean-field band structures are examined. 

PACS numbers: 75.10. Jm,75.50.Ee 



Introduction: The question of the quantum ground 
state of spin- 1/2 Heisenberg Hamiltonian 

H = J2^-^ (1) 

(ij) 

on the pyrochlore lattice has prompted active research 
for nearly two decades 1-6 . A popular scheme, employed 
in Refs. 1-4, was to solve the isolated tetrahedron prob- 
lem exactly and to couple the nearby disjoint tetrahcdra 
in the weak exchange energy J'. A non-magnetic state 
with exponentially short correlation length was identified 
in Ref. 1, and Refs. 2-4 noted the dimer instability in the 
ground state. A similar dimer instability was discovered 
in the large- N approach in Ref. 5. While Refs. 2-4 begin 
with the dimer basis of a single tetrahedron solution to 
carry out perturbation in J', Ref. 5 starts with a trans- 
lationally invariant solution and finds that dimerization 
occurs as a spontaneous symmetry breaking. 

In this paper, we adopt the fermionic mean-field fol- 
lowed by variational Monte Carlo (VMC) treatment to 
address this issue. Prior work in the same spirit for the 
Kagome lattice can be found in Ref. 7. We find chiral spin 
states 8 with non-zero averages of the three-spin operator 
(Si ■ Sj x Sfc) for the elementary triangular unit formed 
by (ijk) sites to be the likely non-magnetic ground state 
of the Heisenberg model realized on pyrochlore lattice, 
under the assumption of uniform nearest-neighbor bond 
amplitudes | (S, ■ Sj}\. The flux through the triangles are 
found to be ±7r/2 at the mean- field level, but reached a 
smaller value after the Gutzwiller projection was carried 
out. 

Mean-field Theory: In rewriting the spin operator 
as a fermion bilinear and introducing the mean-field vari- 
able Xij — {ftfi)i one arrives at the mean-field Hamil- 
tonian 

Hmy = - E E + E V2), (2) 

i j£i i 

with the Lagrange multiplier A,; enforcing the occupation 
constraint at each site. In doing the mean-field calcula- 
tion we drop the spin index and work with the half-filled 



case {ft fi) = 1/2- The summation over all nearest- 
neighbor sites j with respect to i is indicated by Yljev 
To set up the coordinates, we place the four corners 
of a single "up" tetrahedron (uT) at (0,0,0), (1,0,0), 
(1/2,73/2,0), (1/2,1/2^3,^/273), then displace it by 
integer combinations of e\ — (2,0,0), 62 = (1, \/3, 0), 
and £3 = (1, 1/V3, \/8/3) to generate the pyrochlore lat- 
tice: n\e\ + niei + 71363, rii=integers. Each uT is met 
at four corners by down tetrahedra (dT) and vice versa 
as shown in Fig. 1. The self-consistent mean-field cal- 
culations were run for L x L x L = L 3 lattice with 4L 3 
lattice sites. L refers to the number of uT's along each 
e Q direction. The lattice contains an equal number of uT 
and dT blocks. 

A completely unrestricted minimization of (Hmf) re- 
sulted in the ground state with each site paired into a 
dimer while all dimers are disconnected from one another, 
in accordance with Rokhsar's general observation 10 . The 
extensive degeneracy of the dimer ground state will be 
lifted at higher orders in 1/N in a large- N expansion 
to give rise to a ground state with (possibly) restored 
translational symmetry. Furthermore, the fully dimer- 
ized state carries an energy of —0.375 per site, much 
higher than some of the uniform states we consider in 
this paper. Therefore, we move with the idea that the 
true ground state of Eq. (1) is better captured by the 
uniform amplitude ansatz \x%j\ = Xi Du t with arbi- 
trary phases: Xij = X e y ■ F° r a single tetrahedron, 
such a search yielded solutions where the flux through 
the four triangular faces are all equal to $ = +7r/2 or 
$ = — 7r/2. The flux is defined from the directed product 
e »* = e i\<l>i j +<t> j k+4>ki] as tne three sites of a triangle (ijk) 
are traversed in a counterclockwise manner as viewed 
from outside the tetrahedron. For the lattice problem 
with uniform \xij\j the mean-field ground state is found 
to be the one with staggered chirality: $ = tt/2 for all 
the uT's and $ = —tt/2 for all the dT's, or vice versa. 
No solutions were found where different faces of a given 
tetrahedron, either up or down, carried different amounts 
of flux, or flux other than ±7r/2. The flux through the 
hexagons of the pyrochlore lattice were all zero. The 
mean-field ground state solution is in complete confor- 
mity with the "Rokhsar rule" 11 of the flux for different 
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(a) [|, 0]-flux state 




(b) [|,-|,0]-flux state 



(c) [0,0,jr]-flux state 



FIG. 1: (color online) \ij bond configurations producing (a) [§■,■§, 0], (b) , — ^, 0] and (c) [0, 0, n] flux states. Setting | = 1, 
the direction of the arrows indicates \ij = +* going from i to j. The shaded bonds are lying on the back sides. Independent 
sites are labeled 1 through 4 (one up tetrahedron in a unit cell) in (a) and 1 through 16 (four up tetrahedra in a unit cell) in 
(b) and (c). In (c), the thick(dotted) lines represent bonds with \ij = 1( — !)■ 



types of polygons. 

Variational Monte Carlo Calculation: Variational 
Monte Carlo (VMC) calculations of the energies of sev- 
eral mean-field ansatz states were carried out including 
the one found in the mean-field calculation. The states 
are labeled by [uT, dT, H], where the three numbers uT, 
dT, and H refer to the flux through the triangles of the 
up tetrahedra, down tetrahedra, and the hexagons, re- 
spectively. We examine four such states, (i) [0,0,0], (ii) 
[1,1,0], (iii) [f,-f,0], (iv) [0,0,7r]. The mean field re- 
sult corresponds to case (iii). For future reference, we 
denote (ii) and (iii) as uniform and staggered flux states. 
The other two states, (i) and (iv), do not break time 
reversal symmetry. The choice of the trial states are mo- 
tivated by and parallels those in Ref. 7 for the Kagome 
lattice. The Xij bond patterns that generate each state 
at the mean- field level are shown in Fig. 1. With Xij' s 
given in Fig. 1 as input to Eq. (2) one can diagonalize 
.Hmf to obtain the mean-field ground state |'0mf) as a 
Slater determinant. The evaluation of the energy and 
other operators X are carried out in the projected space 
\s) Qs) spans all the states with one spin per site) by 



J 



'MF/ 



,(3) 



D s (V'mf|s)(s|V'mf) (s\iPmf) 



where P(s) = |(s|?/W| 2 / Yls I( s I^mf)| 2 is the probability 
weight used in the Monte Carlo procedure 12 . 

The unit cell includes a single tetrahedron for the 
[0,0,0] and [§, f , 0] states, and 2x2x1 tetrahedra for the 
[§,— f,0] and [0,0, n] states. Periodic boundary condi- 
tions (PBC) generate degenerate states at the Fermi level 
for the [0,0,0] and [§,§,0] cases, which can be lifted 
by applying anti-periodic boundary conditions (aPBC) 
along one of the directions, e.g. e\ , for [0, 0, 0] , and along 
two directions, e.g. e\ and £2, for [f , f ,0]. In the case 
of [■§,— f,0] and [0,0, 7r] the unit cell includes four up 
tetrahedra, two in the ei-direction and two in the 63- 
direction, for the Xij patterns shown in Fig. 1(b) and 
(c). PBC suffices in both these cases since no degeneracy 
occurs at the Fermi level. 



r 



The VMC energies of the trial states we propose are 
listed in Table I. Several independent Monte Carlo simu- 
lations were made with each simulation consisting of over 
10 5 steps (each step means one sweep of the whole lat- 
tice size) for each state to obtain reliable estimates of the 
energy. Due to the rapid increase of the system size with 
L, the calculation is currently limited to L = 4. For L up 
to 4, the two chiral states turn out to have much lower 
energies than the two non-chiral states. The energies of 
the two chiral states are very close. 

Spin-spin Correlation and Chirality: The spin- 
spin correlation function can be computed using the re- 



lation (Si • Sj) = (3/ '4)(crfcr|) due to the spin isotropy 
of the ground state. Table II displays the spin-spin cor- 
relations obtained for all four flux configurations as a 
function of the distance for L = 4. We measure the cor- 
relation along the six directions, ei, £2, £3, ei — 62, e^ — £3, 
£3 — ei, and take averages. The results for any one di- 
rection is consistent with those taken along any other. It 
was difficult to reach a large enough size to discriminate 
an algebraic decay of the correlation against an exponen- 
tial one. The sign of (S^ • Sj) alternated with distance 
as expected in an antiferromagnetically correlated state, 
except for [0, 0, ir] where the signs of the second- and the 
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TABLE I: Energy per site of various states obtained from 
VMC depending on the number of up tetrahedra L in each 
direction. The two chiral states have much lower energies 
than the non-chiral states. Statistical uncertainties lie below 
the digits shown. Boundary conditions, aPBC or PBC, used 
in the calculation are listed. 



L [0,0,0] 


[f.S.o] 


[f,-f,0] [0,0, tt] 


aPBC 


aPBC 


PBC PBC 


2 -0.372 


-0.478 


-0.466 -0.374 


4 -0.374 


-0.459 


-0.456 -0.375 



in Fig. 2. The two states |s+) and |s_) are obtained as 
CCW and CW rotations by one lattice site of the original 
spin configuration |s) with the rotation axis chosen to 
point out of the triangular faces of the tetrahedra. When 
the state \s) contains all three spins up or all down for 
a given triangle, the chirality operation gives zero. The 
average (xm) is obtained from 



<X 



123/ 



(5) 



third-neighbor correlations are reversed (See Table II). 
The fast decay of the spin-spin correlation is consistent 
with an earlier finding of Ref. 1 based on the cluster 
perturbation method. 



Here P(s) is the statistical weight |(s|V')| 2 /EJ( s IV')| 2 
for a given mean-field state \ip) . Following the usual man- 
ner of updating the configuration by Monte Carlo meth- 
ods, one has to calculate the ratio ((s_ \ip) — (s + \ip))/ (s\ip) 
for each state \s). 



TABLE II: The spin-spin correlation function (Si ■ Sj) with 
respect to distance for the four flux states. Five independent 
10 MC steps were used for each data point. 

\i-j\ [0,0,0] [%,%,0] [%,-g,0] [0,0, tt] 

1 -0.1249(2) -0.1533(1) -0.1522(1) -0.1250(1) 

2 +0.0088(3) +0.0194(2) +0.0122(1) -0.0054(3) 

3 -0.0027(1) -0.0129(1) -0.0024(3) +0.0012(2) 

4 +0.0007(2) +0.0060(2) +0.0006(3) +0.0001(6) 

Both the mean-field search and the VMC calculation 
suggest that the chiral states with spontaneously broken 
time-reversal symmetry (T) is vital in the understanding 
of the ground state correlation of the Heisenberg model 
on the pyrochlore lattice. We develop below an extension 
of the VMC method which allows the calculation of the 
chirality and evaluate it for the chiral states. 




FIG. 2: (color online) The two states, |s+) and \s- }, pro- 
duced by the chirality operation on a state \s). 

The scalar chirality X123 = (Si • S2 x S3) is equivalent 
to 8 



1 T - ~ u (X12X23X31 - X13X32X21) 

xlm X12X23X31) = TT- • 4 

2 li 

with Xij = 2~2a fiafjcr- ^ i s tne difference of the cyclic 
permutation of the spins done in the counterclockwise 
(CCW) and clockwise (CW) directions. Unlike in the 
evaluation of the spin-spin correlation where the action 
of erf ct| on a given basis state \s) is diagonal, the outcome 
of Xi23|s) is not proportional to \s) itself. Instead, we 
obtain the relation X123 \s) = i(|s+ ) — |s_ )), as depicted 



TABLE III: The real and imaginary parts of the averages 
of (s+\ip) /{s\tp) and the flux (in units of it/2) through the 
triangle for [?, 5,0] (left three columns) and [■§,— 5,0] (right 
three columns) 



L 


7?123 |Xl23 1 $1231 


??123 X123 |*123 1 


2 


0.376 0.356 0.483 


0.368 0.325 0.461 


4 


0.314 0.358 0.542 


0.312 0.358 0.544 



In the VMC calculation we take averages of 
(s-\ip)/{s\ip) and (s + \i{j) / (s\ip) separately. The two 
quantities turn out to be complex conjugates with ex- 
tremely high accuracy, so one can denote the averaged 
( S+ |V)/<s|V> and (s_|V)/<s|V) as (r7 123 )/2 ± i<Xi23>/2, 
respectively. The flux $123 is deduced from tan(<j>i 2 3) = 
(X123)/ (7)123) • The signs of the flux after the projection 
turned out to be in perfect accord with the mean-field 
predictions. On the other hand, the amount of flux is re- 
duced from the mean-field value it/2 after the projection 
(See Table III). The chirality was zero within a statis- 
tical error for the non-chiral states, [0,0,0] and [0,0, tt]. 
The increase of the average flux <!> with the system size 
is consistent with the scenario of a long-range ordering 
of the chirality in the ground state. 

Band Structure: The band structures of the four 
flux configurations in Fig. 1 have been analyzed along 
the three orthogonal directions k x ,k y , k z as well as along 
k a defined to lie along the three e Q directions, a = 1, 2, 3. 
The relations between the two sets of momenta are k\ = 
k x , k2 — k x /2+^/3k y /2, and k$ — fc a; /2+v / 3fcy/6+v / 6fc z /3. 
Flat bands lying exactly at Ep were observed in the 
[0, 0, 0] and [f , f , 0] states, and above Ep for the [0, 0, tt] 
state. No flat bands exist for [f , — f ,0]. We describe the 
respective band structures in more detail, (i) [0,0,0]: A 
two-fold degenerate flat band lies exactly at Ep for each 
k a direction. The other two non-degenerate bands lying 
below Ep are dispersive, (ii) [§,|,0]: A two- fold de- 
generate flat band at Ep = was observed along the k z 
direction, k x = k y = 0. The other two bands are given 
by ~ ± cos(^/2/3fc 2 ). Along each of the k a directions, 



4 



the uppermost and the lowermost bands (which are re- 
lated by particle- hole symmetry) are flat, and the other 
two dispersive bands cross at Ep as k a equals a multiple 
of 7r. There is no gap for this, or the [0,0,0] flux con- 
figuration. This explains the existence of a Fermi level 
degeneracy in the finite-size mean-field energy spectra, 
(iii) [§,— §)0]: Four doubly-degenerate bands lie below 
the Fermi level and the others lie above it, separated 
by an energy gap. Dispersion along the k a directions 
are displayed in Fig. 3. (iv) [0,0, n]: Of the 16 bands, 
the uppermost band lying above Ep is 8-fold degenerate 
and flat along each of the k a directions. The remaining 
ones are the four, 2-fold degenerate bands lying below 
the Fermi level, and separated by a gap from the upper- 
most one. A gap separates the occupied from the empty 
bands for the [§,— f,0] and [0,0, tt] flux states, which 
also explains the absence of a Fermi level degeneracy in 
an earlier mean-field calculation. 




kj/71 



k 3 /7i 



FIG. 3: Mean-field energy bands for the [f , — f, 0] state along 
k a — k-e a directions, < k a < 2ir. Each band is two-fold de- 
generate. The occupied bands shown here are symmetric with 
the upper, unoccupied bands due to particle-hole symmetry. 
The horizontal line Ef — is the Fermi level. Dispersion 
along other directions (not shown) also show a full gap. 

The presence of flat bands in the [0,0,0] and [§, f ,0] 
flux states over a linear segment of the Brillouin zone 
suggests that an instability might play a role. In the case 
of the uniform flux state, the degeneracy-lifting terms 
are given by the modulation of the bond amplitudes for 
triangles lying in the plane spanned by ei and ei vectors. 
Specifically it corresponds to i\ — > i(x+X) for the (123) 
triangle of the up tetrahedron, and i\ — > i(x~ ^) f° r 
the (123) corners of the down tetrahedron shown in Fig. 
1 (a). The flux through the triangles remain fixed at 



7r/2. We observed that the introduction of nonzero A did 
not decrease the variational energy. Rather there was 
an increase in the third significant digit of the energy 
value for X/ X up to 0.1 and by 2.7% at \/\ = 0.2. The 
slow dependence of the overall energy is partly due to the 
boundary conditions used, which already opened the gap 
in the finite-size single-particle spectra. Another reason 
for the lack of dependence may be that the degeneracy is 
confined to a linear segment of the Brillouin zone, so that 
not enough states are affected by the degeneracy-lifting 
mechanism. On the other hand, the gapful nature of the 
band structures for [§,— f ,0] will guarantee stability of 
the Gutzwillcr-projected state against small fluctuations. 
At the moment we believe both types of flux states to 
have a chance to represent the uniform ground states of 
Eq. (1). 

Discussion: We find that the combined search using 
the fermionic mean-field theory and Gutzwiller projec- 
tion yields chiral spin liquid states (both uniform and 
staggered flux types) with ordered chiralities as the likely 
ground states of the S = 1/2 Hciscnberg spin Hamilto- 
nian on the pyrochlore lattice. Previous theories based on 
the perturbative expansion around the single-tetrahedron 
solution 1-5 did not find such chiral spin states. Given the 
past claims of dimer instability, a simultaneous search for 
a pronounced chirality and dimcr correlations in the ex- 
act diagonalization of Eq. (1) will be valuable in sorting 
out the contending perspectives. The Fermi level degen- 
eracy one finds in the uniform flux state may have inter- 
esting consequences that will need to be explored more 
carefully in the future. 

Note added: After the submission of this article, there 
appeared a preprint 9 (arXiv:0809.0528vl) pertaining to 
the chiral spin state on pyrochlore lattice. Their conclu- 
sions overlap with ours. 
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